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Abstract 

Static, spherically symmetric solutions with regular origin are in- 
vestigated of the Einstein- Yang-Mills theory with a negative cos- 
mological constant A. A combination of numerical and analytical 
methods leads to a clear picture of the 'moduli space' of the solu- 
tions. Some issues discussed in the existing literature on the sub- 
ject are reconsidered and clarified. In particular the stability of 
the asymptotically AdS solutions is studied. Like for the Bartnik- 
McKinnon (BK) solutions obtained for A = there are two differ- 
ent types of instabilities - 'topological' and 'gravitational'. Regions 
with any number of these instabilities are identified in the mod- 
uli space. While for BK solutions there is always a non-vanishing 
equal number of instabilities of both types, this degeneracy is lifted 
and there exist stable solutions, genuine sphalerons with exactly 
one unstable mode and so on. The boundaries of these regions are 
determined. 

1 Introduction 

The last decade has witnessed an intensive study of soliton and black hole solu- 
tions of the Einstein- Yang- Mills (EYM) theory leading to many unexpected re- 
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suits (see ^ for a review) . In particular it turned out that the strong uniqueness 
theorems for such solutions valid for gravitating Abelian gauge fields (Einstein- 
Maxwell theory) cease to hold in the non- Abelian case. Also the stability prop- 
erties of spherically symmetric solutions change, partly due to the breakdown 
of the Birkhoff theorem in the non- Abelian theory. 

Among the many possible extensions of the EYM system studied in the 
literature the simplest one is obtained by adding a cosmological constant A ^Ej. 
As to be expected this changes the global behavior of the regular solutions in 
a characteristic way from asymptotically Minkowskian to that of deSitter (dS) 
resp. Anti-deSitter space (AdS) depending on the sign of A. This change is also 
accompanied by a change of the moduli space of globally regular ('soliton') and 
black hole solutions. While for a given A > the generic solution is singular 
and there is only a discrete set of soliton solutions (characterized by the node- 
number of the YM potential W) there are two types of generic solutions for 
A < 0, the singular ones and asymptotically AdS solutions with continuously 
varying values Woo of Vl^ at r = oo. This change of the topology of the moduli 
space going from negative to positive values of A leads to, what was called a 
'fractal' structure in [3]. 

A further important difference between the dS and AdS case concerns the 
stability properties of the solutions. While it was found that all solutions with 
A > are still unstable like the Bartnik-McKinnon (BK) solutions |S] obtained 
for A = there exist stable solutions with AdS asymptotics for negative cos- 
mological constant [^El. In fact, there seems to be general agreement in the 
literature |3 El 13 El E] that all solutions with a nodeless YM potential W are 
stable. For soliton solutions the question was studied numerically in j4j and 
stability was found for certain solutions with no or one node. An analytical 
proof was given in for black holes with Wh > 1/ \/3 and in lllj for node- 
less solitons, however in both cases assuming A <C 0. How the authors of I^E] 
arrive at their general statement that all nodeless solutions are stable remains 
unclear and in fact we shall present convincing arguments that it is wrong. In 
addition we give a short discussion of the concept of stability employed in the 
above considerations. 

As in the case A > there are two different types of spherically symmetric 
unstable modes, those of 'topological' type (odd parity) and 'gravitational' ones 
(even parity). While for the BK solutions there is always an equal number of 
unstable modes of both types directly given by the number of nodes of the 
solution this is not true for A < 0. 

Our investigation of the stability properties of the soliton solutions is based 
on a general study of the moduli space of the solutions as a function of A. 
In particular we find two discrete sets of lines in the (A, b) plane, (where b 
parametrizes the solutions with a regular origin), separating regions with a 
different number of unstable modes. First there are the lines with vanishing 
Woo for solutions with n nodes. Crossing these lines the number of odd parity 
unstable modes changes by one El- Furthermore there is a discrete set 
of lines of zero-modes in the (A, b) plane separating regions with a different 
number of even parity unstable modes. Besides the region in the moduli space 
without unstable modes there are regions with exactly one unstable mode, the 
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condition for a genuine sphaleron like the traditional one in flat space jl4j . This 
may have some interesting consequences for the application. 

Combining numerical with analytical tools we give a detailed description 
of the boundary between the domain of asymptotically AdS and the singular 
solutions yielding a simple explanation of the 'fractal' structure observed in |^. 
In particular, based on earlier work jl5j on the BK solutions we will derive a 
scaling law for the small A part of the boundary between these two domains. 

For reasons of simplicity we have restricted ourselves to electrically neutral 
soliton solutions, but the generalization to charged solutions and black holes 
should be straightforward. 



2 Field Equations 

The action of EYM theory with the cosmological constant A has the form 

Seym. = ^ / (-^(i? + 2A) - ^F^..F^^^)^gd'x , (1) 

where G is Newton's constant and g is the gauge coupling. In what follows we 
put G = g = 1 choosing suitable units. 

We will be interested in spherically symmetric solutions with regular origin. 
For a spherically symmetric space-time the line element decomposes into 

ds^ = dsl - r^dn'^ , (2) 

where ds^ is the metric on the 2-d orbit space factorizing out the action of the 
rotation group and dQ'^ the invariant line element of S'^. The 2-d metric can 
always be brought to the diagonal form 

dsl = [e^''df-e^^dR^) , (3) 

where t is a time-like coordinate. Here we are only considering time-independent 
solutions and we naturally choose t to be the Killing time. This means that we 
can perform a further dimensional reduction to a 1-d theory. Then the metric 
functions v, A and r depend only on the radial coordinate R. The latter can be 
chosen arbitrarily exhibiting A as a gauge degree of freedom. 

For the SU(2) Yang- Mills field VF^ we use the standard minimal spherically 
symmetric (purely 'magnetic') ansatz 

Wl^^Tadx^" = W{R) {Tide + T2 sin Odip) + T3 cos Odip , (4) 

where Ta denote the generators of SU(2). The reduced EYM action including 
a negative cosmological constant A can be expressed as 



S = -JdR 6^"+^^ (1 + e-'\{r' f + u'{r' 



^ — ^2 + 



(5) 
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Introducing 



N = e~^r' , Hi = re~^u' + N , U = e'^W , (6) 

we write the corresponding Euler-Lagrange equations in first order form 

re-^N' = - iv) iV - , (7a) 

re'^n' = l + 2\A\r^ + 2U^ - , (7b) 



re^^U' = — - + \^N-KjU, (7c) 

2kN - N^ = 2U^ + 1-^ _^L^ + |A|r.2. (7d) 

Tlie last of these equations is the only remaining diffeomorphism constraint. 

We still have to fix the gauge, i.e. choose a radial coordinate. A simple 
choice is to use the geometrical radius r for R (Schwarzschild coords.), yielding 
e"^ = N. Putting fi = N'^, A = e^/N, and A = 0, Eqs. © become equal to 
Eqs. (6) in ,15 . This coordinate choice has the disadvantage that the equations 
become singular at stationary points of r. This problem is avoided using the 
gauge = r introduced in jl5j . which will also be used here. In order to stress 
the dynamical systems character of Eqs. (|6I7|) we denote the corresponding 
radial coordinate by r and derivatives by a dot. We thus get 

r = rN , (8a) 

= K-N , (8b) 

W = rU , (8c) 

N = {k- N)N -2U^ , (8d) 

k = l + 2\A\r'^ + 2U'^ - k'^ , (8e) 

Using the constraint Eq. ((31) the equation for can also be written as 

A = i(l + lAIr^ - ^^^^)) - '-N^ - . (9) 

For A = these equations coincide with the Eqs. (49,50) of JH]- An im- 
portant quantity is the 'mass function' m = r(l + |A|r^/3 — A^)/2 obeying the 
equation 

^=iu'+ ^^l~^^' yN. (10) 

The value of m at r = oo is the AdS mass of the solution. Due to the 
presence of the cosmological constant the vacuum solution obtained for W"^ = 1 
is no longer Minkowski space but Anti-deSitter (AdS) space 

ds^ = ^dt^---r'dn\ (11) 
A' 
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with 

^^{r) = 1 + ^ , (12) 

while the counterpart of the Schwarzschild solution for A < is the Schwarzschild- 
AdS solution obtained for 

M(r) = 1 - — + ^ , (13) 

describing a black hole with AdS asymptotics for Af > 0. 

Likewise there is the Reissner-Nordstr0m-AdS (RN) solution for W = and 

Kr) = 1 + - + 4- , (14) 

r 3 

carrying an Abelian magnetic charge. The function /i has either two simple 
zeros, a double zero or no zero at all for positive r. The double zero occurs at 
Tq = (v^l + 4|A| — 1)/2|A| for a certain value Mo(|A|). As long as fi has zeros 
the solution describes a charged black hole with AdS asymptotics. The zero of 
fj, at the smaller value of r corresponds to the Cauchy horizon of the Reissner- 
Nordstr0m black hole. A double zero of fx is the position of an extremal horizon. 
Solutions without zero have a naked singularity at the origin. 

3 Global behavior 

We want to study the global behavior of solutions with a regular origin. As in 
the case without cosmological constant the field equations are singular for 
r = 0, but using the methods of jl5j it is straightforward to show that there is 
a 2-parameter family of regular solutions with the behavior 

W{r) = l-br^ + Y^H^b - 8b^ + 2|A|)r^ + 0{r^) , (15a) 

^,{r) = 1 - (46^ - + ^b\3b + \A\y + 0{r^) , (15b) 

3 15 

A(r) =^o(l + 46V - —(3&- 186^ + 2|A|)r^ + 0(r^)) , (15c) 

where b and Aq are free parameters. A trivial rescaling of the time coordinate t 
allows to put Aq = 1 and we will assume this condition henceforth for solutions 
with a regular origin. 

We may start with a solution with this behavior at r = and integrate 
the field equations for increasing r. This way we can extend the solution 
until we hit a singular point for some finite value Tg or otherwise to arbitrarily 
large values of r. As long as N{t) is positive r(r) will grow monotonously. 
However, N{t) may have a zero for some value tq. If the function A would stay 
bounded at tq this would lead to a zero of the metric coefficient e'^'^ implying 
a (cosmological) horizon. However, it is easy to see from Eq. (jS^) that this is 
impossible, because k > 1. Thus the function fj, cannot become negative as 
is e.g. claimed in |^. In fact it is only the function that becomes negative 
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and consequently the function r(r) has a maximum and decreases for larger 
r. It follows from Eq. ((HJl) that N, once negative, cannot become positive 
again. As in the case without cosmological constant one can show that tends 
to — oo for some finite value Tg and r runs back to zero. There the solution 
becomes singular in complete analogy to the case A = |15j . since the terms 
proportional to A may be neglected for r — > 0. One finds a Reissner-Nordstr0m 
type singularity of the metric [Ej 



W(r) = Wo + ^jYT^r' + War' + O(r') , (16a) 



,M = M^-2^ + 0(I), (16b) 

A{r) = Ao + 0{r^) . (16c) 

with the arbitrary parameters Wq, W3, Mq and ^o- Topologically the 2-spheres 
of fixed r of these solutions foliate a 3-sphere with a singularity at its second 
'origin', suggesting to call them 'Singular Compact' solutions' (SC). From the 
number of parameters it follows that this is a 'generic' class of solutions. 

The situation where N{t) stays positive is a bit more involved and we 
will have to distinguish three cases. Since r(r) is monotonously growing it 
will either have a finite limit or it will tend to infinity. In the first case the 
situation is practically the same as without cosmological constant. Following 
the arguments given in [H] one concludes that the solution exists for all r and 
tends for r — > 00 to the fixed point corresponding to the degenerate horizon 
of the extremal RN black hole with (r, W, U, N, k) = (ro, 0, 0, 0, kq) where Tq = 
(v/l + 4|A| -1)/2|A| and = ^1 + 4|A|. Linearizing the Eqs. around this 
point exhibits the hyperbolic character of the fixed point with eigenvalues kq, 
-2, -k;o(1±7)/2 where 7^ = 1-4/^1 + 4|A|. As long as |A| < 15/4 there are 
two complex conjugate eigenvalues leading to damped oscillating eigenfunctions, 
while for |A| > 15/4 these eigenvalues are real and negative. As in the A = 
case the space-time splits into two separate parts. The exterior part is the 
exterior (to the degenerate horizon) of the extremal RN-AdS black hole and 
is geodesically incomplete, while the interior part containing the origin is a 
geodesically complete manifold with the topology of AdS and not a space with 
a boundary as claimed in |7j. 

In the second case r — > 00 there remain two possibilities; either the solution 
is asymptotically AdS with finite mass M = limr^oo ?7i(r) or it is a new singular 
solution, with the mass function m{r) diverging as r^. The AdS type solutions 
have the asymptotic behavior (we use again Schwarzschild coordinates) 

m0^^o..^ ^"-^^-^-^^ l.o(J,), (ira) 

Kr) = l- — +^-^r' + iiW^'-lf + l\Mc')^,+0(^) , (17b) 



Mr) = A^l- — +0 -,]], (17c) 
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with the arbitrary parameters Woo, c, M, and A^o- Like the 'SC solutions this 
is a generic class of solutions. From the linear growth of = -^//I with r it 
follows that T remains finite for r — > oo. 

Apart from this regular AdS behavior at infinity there is a singular class of 
solutions (called SAdS) with W, U, N, and k growing linearly with r. For their 
description we introduce new variables 

W = ^ , fi = -^ , and V = JiW . (18) 
In terms of these variables the field equations become 

rW'='^-W, (19a) 

1 2V^ 
rV' = W(W^ - ^) -2V - , (19b) 

OT/2 _ 1 1 

rfl' = -3fl-^-{W^-^f + \A\ + ^. (19c) 

The singular solution corresponds to a hyperbolic fixed point of these equa- 
tions for r — > oo with finite values Wo,Vo,flQ determined by the vanishing of 
the r.h.s. These values can be expressed in terms oi x = Wq obeying the cubic 
equation 

+ + 3x - 2| A| (1 + x) = , (20) 

which has always one positive root. The values Vq and /zq at the fixed point 
can be expressed in terms of x as Vo = ix3/V(l + x) and /2o = ^x/{l + x). 
Introducing w = W — y/x, v = V — Vq, and x = A — Ao and linearizing Eqs. (|19j) 
in these shifted variables we obtain 

, 2(1 + x) 2(1 + x) , , 

rw' = -w + — -V- ^ ^ 21a 

X Jx 



r't;' = 3tt; - 2(1 + 3x)t; + 4xyxx , (21b) 
rx' = — 4x-v/xti; — 4v^t; + (2x — 3)x • (21c) 

Studying the characteristic equation of this linear system one finds that its 
eigenvalues are always real for x > 0. Since for x — > the eigenvalues become 
1, —3, —4 and the determinant 16x^ + 40x2 + 32x + 12 is positive for x > there 
is exactly one unstable mode for any |A| > 0. 

Our numerical analysis shows that any solution with regular origin belongs 
to one of the types described above and it would be desirable to prove this 
analytically. All the considered solutions have some well-defined asymptotics 
and thus there is no indication of any chaotic behavior. 



4 Moduli space topology 

Next we want to describe the moduli space of the solutions with regular ori- 
gin as obtained from a thorough numerical analysis supported by analytical 
arguments. Our numerical results are obtained using a shooting technique 
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starting from r = (not from e.g. r = 10~^) combined with integral equations 
if necessary as described in |16j . For given A the solutions depend only on the 
parameter b. In contrast to the case without cosmological constant there are 
now two different types of 'generic' solutions, the 'Singular Compact' ones run- 
ning back to r = already present for A = and those with AdS asymptotics. 
Therefore we expect certain open sets in the (A, b) plane corresponding to these 
two types. Fig. Ogives a schematic description of these open sets. The domain 
called 'AdS' corresponds to the AdS type solutions. It has a rather complicated 
boundary, consisting of the curves 'RN' and 'SAdS'. The latter consists of the 
lower boundary and infinitely many arcs connecting the points (A, b) = (0, 
where the 6„'s are those of the n*^ BK solution.^ The set 'AdS' may be further 
subdivided according to the number of nodes of W as indicated in Fig. ^ 

The complement of the closure of 'AdS' corresponds to the 'SC solutions. 
It consists of the disconnected pieces ^Sn {n = 0,1...), containing solutions 
with n zeros of W and iWol > 1- Furthermore there is the set containing 
solutions with any number of zeros but with \Wq\ < 1. These sets are in close 
correspondence with those introduced in ^Sj. Since the intervals between the 
bnS on the A-axis belong to the sets '5„' and the latter are open, they extend 
into the region A 7^ 0, but as will be explained later, their boundaries become 
numerically very quickly indistinguishable from the A-axis for larger values of n. 
Varying the parameter b for fixed A closer and closer to the axis one intersects 
more and more of the sets ^Sn alternating with intervals in 'AdS'. This explains 
the 'fractal' structure close to the A-axis observed in ^1. 

The left and lower boundary ('SAdS') of 'AdS' corresponds to the singular 
solutions with unbounded W{r). Hence Woo in 'AdS' resp. Wq in 'S'n' becomes 
arbitrarily large approaching these boundaries. 

The upper boundary curve 'RN' of 'AdS' corresponds to solutions running 
into the RN fixed point described in the previous section. It starts at the 
value boo ~ 0.7642 of yielding the limiting solution of the BK family for 
A = 0. For |A| < 15/4 the corresponding solutions have infinitely many zeros 
and hence the neighborhood of this piece of the boundary contains solutions 
with arbitrarily many zeros. 

Figs. IS and 13 show the same moduli space as Fig.^ but this time with actual 
numerical data. The dotted curves describe solutions with constant values Woo 
resp. Wq. All the Woo curves emanate from the BK values (0, 6„) on the A-axis. 
Crossing the dashed-dotted lines Woo = the number of zeros of W changes 
by one. The lowest of these lines separates solutions without zero from those 
with exactly one. It bifurcates with the 'RN' boundary curve at the point Pq 
with |A| ~ 5.0646, while all the other zero lines bifurcate with 'RN' at Pqo with 
|A| = 15/4. 

It is easy to understand that the region where b ^ belongs to '^oo'. 

There the A term can be neglected and the situation is the same as for A = 0. 
The same holds for b <C — \/|A| corresponding to the region 6 <C for A = 0. 
In the next section we shall derive more detailed properties of the boundary 
curves. 

^We include the trivial solution W = 1 obtained for feo ~ 0. 
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5 Boundary Curves 



The boundary curves are given by non-generic solutions corresponding to fixed 
points with one unstable mode. Hence they have to be determined by fine- 
tuning one of the parameters b or A. 

In order to determine the 'RN' curve one takes solutions running into the 
RN fixed point. As mentioned in section |31 such solutions are oscillating as long 
as |A| < 15/4. For |A| > 15/4 these solutions can have only a finite number of 
zeros. In fact, our numerical analysis shows that they have exactly one zero for 
15/4 < |A| < Li with Li 5.0646 (the point Pq of Figs. [Hand and no zeros 
for larger values of |A|. Therefore all the curves Woo = in 'AdS' with more 
than one zero merge with 'RN' exactly at |A| = 15/4 (the point Poo), while the 
lowest curve with one zero merges with 'RN' at Pq. The curves for nodeless 
solutions with < Woo < 1 run out to |A| = oo eventually approaching the 
boundary 'RN'. 

Next we shall consider the asymptotic behavior of the boundary for large 
values of |A|. As a first step we rescale r — > f = |A|^/^r in order to remove 
the explicit A dependence from the field equations. To retain the expansions 
Eq. (|15b.c) at r = we have to compensate this by a corresponding rescaling 
of 6 — 6 = |A|^^/^6. In fact, the expansion Eq. I!\15h ) for W even requires to 
put b = 1/2. Thus asymptotically the boundary of 'AdS' becomes the parabola 
|A| = 46^. Fig. |1] shows the remarkably perfect fit of the boundary with the 
(empirically adapted) parabola |A| = 46^ — 3b — 0.1 even for small values of |A|. 

The asymptotic form of the solution on 'RN' can also be obtained rather 
simply. Introducing W = (W — l)/r we find 

W = -r/2 + 0{\A\-^/^) , (22a) 
U = -m + 0{\A\-^/^) . (22b) 

Neglecting non-leading terms the Eq. @ for N becomes 

N = fN^ = -(1 - iV^ - 2[/2 - AW^ + f2) = -{1 - - 2r'^N'^) . (23) 

This linear equation for N'^ = fi can be solved imposing the boundary condition 
/i(0) = 1 with the simple result 

H=- / e" dr' . (24) 

r Jo 

For f — > cxD the solution behaves as = f^^/\/^+ 0{f^^) and U = —1/^/2 + 
0{f-^). 

It seems to be much more difficult to obtain an asymptotic form of the SAdS 
solutions for |A| oo as compared to the solutions running into the RN fixed 
point. 

The left boundary of 'AdS' with its infinitely many arcs looks at first sight 
rather complicated. However, making use of the results of ^1,5, we shall derive 
a simple scaling law for the asymptotic form of these arcs for large n.^ The 

similar argument was used in [H| to derive a scaling law for the ADM mass of the 
solutions. 
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basic idea is to relate the neighborhood of the points (0, 5„) to that of the point 
(0, bo = 0). As a first step we consider the solutions with |A| <C 1 and \b\ <^ 1. 

A suitable rescaling in this case is again f = \/|A[r. Putting U = |A|~^/^C/ 
and keeping only leading terms in |A| we obtain from Eqs. (jSJ 

r = fN , (25a) 

N = {k- N)N , (25b) 

K = 1 + 2f2 - , (25c) 

W = rU, (25d) 

U = W{W^ -l)/r + {N - k)U , (25e) 

together with the constraint 

2kN - iV^ = 1 + f2 . (26) 

Since the equations for N and k, are free of matter terms we obtain the AdS 
metric. The remaining Eqs. ()25tl.e) for W and U are just those for the YM field 
in the AdS background with |A| = 1. Clearly we have to rescale also b = \A\b 
to obtain the expansion W = 1 — bf^ + 0{f^) at the origin. 

Numerically one finds a finite interval bo <b <bi with bo ~ —0.294582 and 
6i PS 2.970497 for which W stays bounded. Outside this interval W diverges for 
finite f. There is an exact solution UH] W{f) = (l + f^/3)~^/^ corresponding to 
b = 1/6, which separates solutions with and without zero, solutions for 6 < 1/6 
having no zero, while solutions with b > 1/6 having exactly one. 

Next we analyze the boundary in the neighborhood of the BK points (A, b) = 
(0, bn)- In 15 it was argued that for b ^ bn there is some r„ ~ roe such that 
for r rn the solution can be parametrized in the form 

N{r) = l-^^^^, (27a) 
r 

W{r) = {-ir{l-bt\r,by-'-^) , (27b) 

with slowly varying functions M„, b^\ and c^. While the coefficients of the 
convergent modes Mn{b,r) and Cn{bn,r) have non-zero limits for r ^ oo that 
of the divergent mode bl^\bn,r) vanishes in this limit. However, assuming 
transversality its derivative dbi^\bn,r)/db has a non-zero limit b^'^ and thus 
6^°°^ fa b^'^{b — bn)- Taking r ^ Vn and b — bn <C 1 we can neglect the 
terms with M„ and Cn still keeping <C 1. This implies the relations 

~ 1 and U/2 « (W it l)/r characterizing the stable manifold of the singular 
point r = 0, i.e. valid near r = for solutions with a regular origin. Thus 
we succeeded to relate solutions with 6 ~ 6„ to solutions with 6 ~ 0. All this 
remains true for non-zero |A| <^ 1 and we can repeat the arguments from above 
for small A and b. Rescaling again r f = |A|^/^r we will find the same 
behavior of the solution with b ^ bn for r ^ r„ as with 6 ~ for f > if 
b = 6i°°^|A|~^ — 6„)|A|^-'^. Arguments analogous to the ones used in 

[15j to derive scaling laws for the parameters of BK solutions with large n show 
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2mr 



that b^'^ scales like ~ const. • e . Therefore we expect to find the same 
linear relation b = 6o,i|^| valid for the boundary curves near {A,b) = (0,0) 

2n7V 

also near the points (0, if we rescale A by the factor e ~ (600 — b) . 
Fig. [21 shows this scaling law is in fact well satisfied already from n = 3 on 
and even holds not only for the whole arcs joining subsequent 6„s, but also the 
neighboring curves of constant Woo resp. suitably rescaled Wq. 

The scaling argument for the singular solutions parametrized by Wq has to 
be different from that for the AdS ones. The singular solutions start from a 
regular origin r = 0, have a maximum of r ('equator') at some point r = and 
then run back to r = 0, where N —00. In ^5] an asymptotic formula for re, 
We, and [/p. for small values of b was derived for the case A = 



1 1 1 2 VF^ 

re ^ f\b\~2 , \We\ ~ 2-633r| , \Ue\ ~ -j^ , (28) 

V2re 

where f ~ 5.317 for b > and f ~ 1.746 for 6 < 0. Observe that We grows as 
for b ^ 0. The derivation of these formulae is based on two steps. First 
the equation for W is integrated on a flat background, leading to a divergence 
of \W\ « 1/(too — r) for some finite value Tqo and a corresponding finite value 
Too = f|6|~^/^. The scaling with is a consequence of the scale invariance 

of the flat YM equation. The two values of f for the two sign choices of b 
are determined by a numerical integration of the flat YM equation. Using this 
diverging solution for W, the equation for is then integrated from A^ = 1 to 
A^ = 0, the equator, yielding flnite values Te < Too, We, Ue, and rg. Since this 
change of A^ happens on a r interval whose length tends to zero for 6 ^ 0, we 
get Te ^ . This derivation remains essentially valid taking into account the 
cosmological term, but there are two modiflcations. Firstly, the function W is 
now considered in the AdS background. This step yields functions f{b) (with 
b = b/\A\) for 6 > resp. b < 0, replacing the two values of f from the flat 
case (compare Fig. |H1 for a numerical determination) . Secondly, the function N 
has now to run between the AdS value ^1 + |A|r^/3 and the equator A^ = 0. 
Hence we obtain 

61 1 1 2 lAIr^ 1 W'^ 

re^r{^^)\bn, |W,|« 2-633r|(l + U^)3 , \Ue\ ^ ^ , (29) 

From a systematic point of view it is better not to stop at the equator as 
was done in but to continue the solution all the way back to r = and 
to characterize it by its value Wq at r = 0. Therefore we decided to plot Wq 
instead of We (compare Fig. [51). For small values of b one finds Wq ^ We, since 
Wo scales like l^l"^/^ for 5 — > (remember We ~ l^l""*^/^). In order to describe 
the growth of W between the equator and the singular origin we rescale all the 
variables W ,U, r, N, k and also the independent variable r by W ^ W|6|~^/^ 
etc. and use approximate equations derived from Eqs. © keeping only leading 
terms for 6 — > 

f = fN, (30a) 
W = ftj , (30b) 
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N = {R- N)N - , (30c) 
R = 2U'^ -r"^ , (30d) 
- - 

U = — + (N-R)U, (30e) 
r 

Observe that there is no more dependence on A in these equations. In the 
hmit 6 — > the rescahng stretches the finite r interval to infinite length and 
the equator becomes a fixed point of the approximate equations. Because of 
the different scaling at the equator all the rescaled variables W etc. vanish at 
this fixed point like negative powers of r with fixed coefficients. The only free 
parameter (apart from the sign of W) is the value re of r. Integrating the 
Eqs. (|3fl|) numerically over the infinite r interval from r = rg to r = we find 
\Wq\ ~ 1.14871re. Using Eq. H29|l we obtain the simple relation 

i_ 

\Wo\ = 1.14871 ^ ll= =\w/^ . (31) 

The argument about the mapping of the neighborhood of the points (0, b>^^ 
to that of (0, 0) by a suitable scaling works as for the asymptotically AdS 
solutions with one difference. In contrast to Woo the quantity Wq has to be 

(c>o) — "^/^ ^ 

rescaled, since Wq ^ bn ~ e ^ . Correspondingly the Wq curves in FigO 

refer to Wq values rescaled by a factor ev^ ~ 6 from one Sn domain to the 
next. Since the relation (|3T|) does not refer to n it should hold equally well 
in all sectors Sn- Since it is derived for \Wo\ » \We\ ^ 1 it is easier to test 
for larger n because of the rescaling. Fig|3 displays how well this relation is 
satisfied for some selected values of Wq and n. 



6 Stability 

While the BK solutions as well as their counterparts for A > studied in 
El Cni are known to be unstable [201 El] it was observed in !11 El Ej that 
solutions in a certain sub-domain of 'AdS' are stable. Actually, as for BK 
solutions there are two different kinds of instabilities. There are instabilities 
of 'topological' type ('odd parity' sector) comparable to those of the flat space 
sphalerons J7T and 'gravitational' ones ('even parity' sector) jl2l l2Uj . As was 
shown in jl^ for the BK solutions and generalized to solutions with A > 
in the number of unstable modes of the topological type is equal to the 
number of zeros of W. Although the same seems to be true empirically in the 
BK case in the even parity sector there is no general proof (although there is 
a rather convincing argument ^2)- these rules would be true also for A < 
one would expect stability in both sectors for solutions in 'AdS' without zero of 
W, as was claimed in 4 . Although the proof for the odd parity perturbations 
can be generalized to A < [3 ^j, the situation for the even ones is more 
subtle. In fact, we shall provide numerical evidence that the claim in [H El El 
is wrong and not all the nodeless solutions are stable. 

Before entering the detailed discussion of our numerical results, we would 
like to discuss some more general aspects of the stability issue. Stability is 
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meant here, as well as in the literature mentioned above, to be stability on the 
linearized level, i.e. under infinitesimal perturbations. In the present context of 
a stationary background the perturbations are conveniently decomposed into 
Fourier modes with respect to the Killing time t and stability means absence of 
imaginary frequencies, i.e. exponentially growing modes. The linear differential 
equations for these Fourier modes contain in general some gauge degrees of 
freedom (independently from the gauge chosen for the background solution), 
which have to be fixed. Clearly the gauge modes have to be separated out and 
the physical spectrum should be gauge independent. However, the structure 
of the resulting linear system may be quite different and favor some particular 
gauge choices. Apart from that there is a choice of boundary conditions to 
be made. The latter has two basically different aspects. On the one hand 
there may be physical conditions for the behavior of the perturbations on the 
boundary, on the other hand one usually requires self-adjointness of the linear 
differential operator. The latter is a technical assumption made in order to 
have an expansion theorem for the general class of perturbations considered. 
In the case of asymptotically fiat solutions one may restrict the consideration 
to perturbations of compact support, although the eigenmodes themselves will 
generally not have this property. In this case it is important to have essential 
self-adjointness of the linear differential operator on the perturbations with 
compact support in order to fix the spectrum. 

Let us now turn to the system under consideration. We shall restrict our- 
selves to the spherically symmetric perturbations of even parity. It is conve- 
nient to use the gauge where the variations are taken at fixed r ('Schwarzschild 
gauge'), but observe that this does not require the use of Schwarzschild coor- 
dinates. In fact, it turns out to be convenient to use a coordinate p, for which 
the 2-d metric is conformally fiat putting dr = Afidp. With these conventions 
the spectral perturbation problem can be reduced as for A = to a single 
'Schrodinger equation' [£j for the perturbation 5W of W 



with the potential 

u{p) = y—^ - 4— (i - - ,2 ) + 8^ ,3 j ' 

(33) 

where W' = dW/dr. 

Next we turn to the boundary conditions. In contrast to the asymptotically 
fiat case, where p runs from to infinity it runs only over a finite interval 
(0,/9m] in the AdS case. While the potential U 7^2/ p^ at p = is singular it is 
easily seen to be finite at pm- Correspondingly we are in the 'limit point' case 
at /O = 0, but in the 'limit cycle' case at pm according to Weyl's classification 
|23j . This means that the differential symbol does not define an essentially 
self-adjoint operator, but there is a 1-parameter family of boundary conditions 
at Pm- This is related to the fact, that there are two arbitrary parameters 
for the background solution at r = cxd (compare Eq. (|17j) ). namely Woo and 




(32) 
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c. The standard choice in the hterature is 6Woo = — a choice that looks 
reasonable, although not much more physical than any condition of the form 
cos a 5Woo +sin a6c = 0. Clearly the spectrum of the linear differential operator 
depends on the boundary condition. Thus the class of perturbations vanishing 
in a neighborhood of Pm do not determine the spectrum uniquely. After all, this 
is perhaps not too surprising for a space-time, in which time-like geodesies reach 
r = oo in a finite coordinate time. Anyhow, in order to be able to compare our 
results to those in the literature, we make the choice 5Woo = 0. 

A standard way of analyzing a change of stability is to search for zero modes 
within a continuous family of solutions jj^]- Prom Fig. [21 it can be seen that all 
the curves with a fixed value of Woo / with at least one node joining the points 
{0,bn) with (0, have turning points with a maximal value of |A|. These 
points naturally correspond to solutions supporting a zero mode. Assuming that 
the number of unstable modes changes only at these points we get a relation 
between the number of unstable modes of the corresponding BK solutions and 
the branches of the Woo curves. Of particular interest are the curves joining 
6o = and bi. As was already mentioned above, these are just the curves 
with Woo < 0, while the Woo = curve bifurcates with the 'RN' boundary 
for |A| « 5.0646 and those with Woo > run out to |A| = oo approximating 
the upper resp. lower boundary for Woo < 1 resp. Woo > 1- At first sight this 
seems to support the hypothesis that the Woo > curves have no turning points, 
implying the stability of these solutions. A careful numerical analysis, however, 
reveals that this is wrong. In fact, as Fig. |H1 shows the Woo = curve still has 
a turning point at |A| ~ 5.104. Thus the solutions on the upper branch of this 
curve between the turning point and the bifurcation point are expected to be 
unstable. A numerical study of the solution of the linearized equations with 
zero eigenvalue supports this result. Fig. El shows two different 'zero-energy 
wave- functions' 5W for solutions with Woo = taken below and above the 
dashed stability line in FigjHl The lower one has no zero, the upper one has 
one, implying stability resp. existence of one negative eigenvalue corresponding 
to an unstable mode j2ll. Neighboring curves with sufficiently small positive 
values of Woo < 0.0005 also have a turning point and follow the Woo = curve 
close to the 'RN' boundary, but then turn around once more and run out to 
large values of |A|. 

Our discussion shows that there is a region in the moduli space where the 
solutions are stable. It lies to the right of the lowest Woo = curve and of the 
curve for the lowest zero mode and is denoted by A in Fig.jHl Then there are two 
regions B resp. C with solutions that have exactly one even resp. odd unstable 
mode, while in D they have one unstable mode of either parity. Similarly there 
are further regions with 2n — 1 resp. 2n unstable modes between the two types 
of curves for solutions with n nodes. All these curves terminate at the point 
Poo on the 'RN' boundary. 
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Figure 2: The numerically determined moduli space: solid and dashed-dotted 
lines as in Fig. ^ zero-mode curves are dashed, curves for several Woo resp. 
Wq =const. are dotted. 
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Figure 3: Small A part of Fig. |2lwitli a logarithmic scale. 
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Figure 4: The parabola fit for large |A| to the boundary curves of AdS. 
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(b„-b)2 

Figure 5: The rescaled regions Sn- the dashed-dotted hnes are curves W^o = 0, 
the dotted hnes in AdS are curves for |Woo| = 10, those in Sq, . . . , S4 are curves 
for |Wo| = 6.5, 40, 245, 1505, and 9231. 
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Figure 6: The radius r^o for singular solutions in the AdS background as a 
function of 6 = ^; the vertical dotted lines at 6 = -0.294582 and 2.970497 
represent the boundary for singular solutions, the horizontal dotted lines are 
the BK values f = 1.74575 and 5.31661 for A = 0. 
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Figure 7: A plot of the quotient Q of th l.h.s. and r.h.s. of Eq. for solutions 
with iV = 0, . . . , 3 nodes and various values of Wq. 
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Figure 8: Enlarged view of the moduli space near Pq with Wqo curves in the 
range ±4 • 10~^; the regions A, . . . , D with different numbers of even and odd 
unstable modes are separated by the dashed-dotted VFqo curve and the dashed 
zero-mode curve. 
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Figure 9: Zero energy wave functions 6W for nodeless solutions with Wqo = 
for b = 1.5, A = -5.075439 lying below and b = 1.553, A = -5.101469 lying 
above the zero-mode curve. 
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